Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 



(MN WF^ style file v2.2) 



Modeling non-linear effects in the redshift space two-point 
correlation function and its implications for the pairwise 
velocity dispersion 



in 
o 
o 



Biswajit Pandey* and Somnath Bharadwajf 

Department of Physics and Meteorology 
and 

Centre for Theoretical Studies 
IIT Kharagpur 
Pin: 721 302 , India 



0^ 



> 
O 

m 
o 

o 
o 



2 February 2008 



ABSTRACT 

The anisotropies in the galaxy two-point correlation function measured from redshift 
surveys exhibits deviations from the predictions of the linear theory of redshift space 
distortion on scales as large as 20 Mpc where we expect linear theory to hold in 
real space. Any attempt at analyzing the anisotropies in the redshift correlation func- 
tion and determining the linear distortion parameter f3 requires these deviations to 
be correctly modeled and taken into account. These deviations are usually attributed 
to galaxy random motions and these are incorporated in the analysis through a phe- 
nomenological model where the linear redshift correlation is convolved with the ran- 
dom pairwise velocity distribution function along the line of sight. We show that a 
substantial part of the deviations arise from non-linear effects in the mapping from real 
to redshift space caused by the coherent flows. Models which incorporate this effect 
provide an equally good fit to N-body results as compared to the phenomenological 
model which has only the effect of random motions. We find that the pairwise velocity 
dispersion predicted by all the models that we have considered are in excess of the 
values determined directly from the N-body simulations. This indicates a shortcom- 
ing in our understanding of the statistical properties of peculiar velocities and their 
relation to redshift distortion. 
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1 INTRODUCTION 

Galaxy redshifts are not perfectly described by pure Hubble's law. Density fluctuations induce peculiar velocities relative to 
the general Hubble expansion. The peculiar velocities perturb galaxy redshifts which in turn affects their inferred distances 
and this leads to a systematic distortion in the clustering pattern of galaxies in redshift space. The peculiar velocities cause 
the two-point correlation function in redshift space 5" (s) to be anisotropic i. e. it depends separately on the component of the 
pair separation s parallel (s||) and perpendicular (sx) to the observer's line of sight fi. There are two characteristic effects of 
peculiar velocities . On large scales structures are compressed along the line of sight due to coherent ffows into over dense 
regions and out of under dense regions, thereby amplifying ^''(s||,sx). On small scales C^(s||,sx) is suppressed due to the 
structures being elongated along the line of sight by random motions in virialized clusters. 

Kaiser (1987) first quantified the correlation anisotropy that results from large-scale peculiar flows in terms of the power 
spectrum of galaxy clustering. Using linear theory and the plane parallel approximation he showed that the power spectrum 
in redshift space -Ps(k) and it's real space counterpart Prik) are related as 
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P.(k) = + P,(fc) (1) 

where Hk is the cosine of the angle between k and the line of sight n, and /3 ~ /h is the linear distortion parameter. Here 
fim is the cosmic mass density parameter and b is the linear bias parameter which differs from unity if the galaxies represent 
a biased sample of the underlying dark matter distribution. It may be noted that the factor relates peculiar velocities to 
density density fluctuations (Peebles 1980). This is slightly modified in the presence of a cosmological constant (Lahav et al. 
1991) and it is more accurate to use (3 = f{Q,m,)/b where /(f2m) = f^m* + ^I-*^ ~ |f^m(l + ^m)] ■ The important point is that 
the anisotropics observed in Ps(k) can be used to determine the value of (5, and thereby place interesting constraints on the 
density parameter Om and the bias b. This has been the single most important motivation for a substantial amount of the 
research which has been carried out in trying to understand and quantify the nature of redshift space distortions. 

Hamilton (1992) translated Kaiser's linear formula from Fourier to real space. He showed that it is most convenient to 
parameterize the anisotropy of ^°(s|[,sx) in terms of spherical harmonics as 

oo 

e{sn,s±) = Y,^i{s)Pi{^,) (2) 



where s = + s^, jj, = sy/s, Pi(At) are the Legendre polynomials and ^i{s) are the different angular moments of the 

redshift space two-point correlation function. Only the first three even angular moments, namely the monopolc £,o{s), the 
quadrupole ^2(5) and the hexadecapole ^4(5) are non-zero and these can be expressed in terms of the real space galaxy 
two-point correlation ^(r) and its moments which are defined as 

Us) = ^ Tmy-'dy. (3) 
* Jo 

through 

^o{s) = (1 + ^(3 + ^(3') as) (4) 
6(s) = (^/3+^/j') [as)-6(s)] (5) 

U{s) = m + ^Hs) - lUs)] (6) 

The linear analysis predicts a negative quadrupole (i.e. ^2(5) < 0) arising from the squashing of large scale structures along 
the line of sight. 

Hamilton proposed that the observed redshift space correlation function be decomposed into spherical harmonics, and 
the ratio 
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1 + 1/3 + i/32 



(7) 



which is expected to have a constant value (shown in [...] in eq. 7) be used to determine the value of /?. Alternatively, if the 
real space correlation function has a power law behaviour ^(r) oc r '', the ratio ^2(s)/^o(s) is also expected to be a constant, 
and this can be used to determine the vahio of B. 

Investigations using N-body simulations to study the redshift space two-point correlation (eg. Suto & Suginohara 1991, 
Fisher et al. 1994, Brainerd et al. 1994, Bromley, Warren & Zurek 1997) find deviations from the linear predictions out to 
scales as large as 20 Mpc and even larger where linear theory is known to be valid in real space. Such deviations arc also 
seen in the redshift space two-point correlations determined from different redshift surveys (eg. Landy, Szalay & Broadhurst 
1998 (LCRS), Peacock et al. 2001 (2dFGRS), Hawkins et al. 2002 (2dFGRS)). In addition to the squashing predicted by the 
linear analysis, the two-point correlation function exhibits an elongation along the line of sight at scales as large as 20h~^ 
Mpc. This causes the quadrupole moment to remain positive even at scales where one would expect linear theory to hold in 
real space. The values of Q{s) which are expected to be constant (eq. 7) do not flatten out to scales as large as 20 Mpc 
in N-body simulations, nor is the flattening observed at these scales in the redshift surveys. All this indicates that there arc 
non-linear effects which are important in the mapping from real space to redshift space at length-scales where linear theory 
is known to be valid in real space. 

The most popular approach is to attribute the deviations from the linear predictions to the effects of the random peculiar 
velocities of galaxies located in virialized clusters and other highly non-linear regions. This effect is incorporated through 
a phenomenological model (eg. Davis & Peebles 1983, Fisher et al. 1994, Peacock & Dodds 1994, Heavens & Taylor 1995, 
Marzke et al. 1995, Ballinger.Peacock & Heavens 1996, Tadros & Efstathiou 1996, Bromley, Warren & Zurek 1997, Ratcliffe 
et al. 1998) which assumes that at large scales the deviations from linear theory can be incorporated by convolving the linear 
redshift space correlation function ^f, with the line of sight component of the random, isotropic pairwise velocity distribution 
function f{v). The resulting non- linear redshift space two-point correlation function is given by 
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r(sii,s±) = J aisii+v,s^)f{v)dv (8) 

where the distribution function /(v) is normahzed to f {v)dv = 1. The authors who have invoked this model have generally 
adopted either a Gaussian or else an exponential pairwise velocity distribution function. In both cases, the distribution function 
has only one unknown quantity aji which is the velocity dispersion of the random component of the pairwise peculiar velocity 
of the galaxies. In this model, the observations of the anisotropies in can be used to jointly determine the value of j3 and 
aR. This has recently been accomplished for the 2dFGRS where they find (3 = 0.49 ±0.09 and an = 506±52kms~^ (Hawkins 
et al. 2002). 

An alternative approach is to attribute the deviations from the linear predictions in to non-linear effects arising from 

the coherent flows. Taylor & Hamilton (1996), Fisher & Nusser (1998) and Hatton & Cole (1998) have used the Zel'dovich 
approximation to analytically study the behaviour of the redshift-space power spectrum in the translinear regime. They find 
that the results from the Zel'dovich approximation are in reasonable agreement with the predictions of N-body simulations, 
indicating that the coherent flows may be making a significant contributions to the non-linear effects observed in the redshift 
space two-point correlation function. 

In a different approach to studying the deviations in from the linear predictions at scales where linear theory is known 
to be valid in real space, Bharadwaj (2001) has considered the non-linear effects introduced by the mapping from real space to 
redshift space. Under the assumption that linear theory is valid in real space and that the density fluctuations are a Gaussian 
random filed, has been calculated taking into account all the non-linear effects that arise due to the mapping fi:om real to 
redshift space. It may be noted that the original calculation of Kaiser (1987) and Hamilton (1992) treats the mapping from 
real to redshift space to linear order only. 

In summary, at large scales where linear theory is known to be valid in real space, the commonly used phenomenological 
model for attributes all the deviations from the linear predictions to the effects of random motions on the mapping from 
real to redshift space. On the other hand, Bharadwaj (2001) calculated incorporating all the non-linear effects which arise 
in the mapping from real to redshift space assuming that they are caused only by the coherent flows. In all probability, the 
deviations from the linear predictions found in ^" in the N-body simulations and in actual redshift surveys is a consequence 
of non-linear effects in the mapping from real space to redshift space arising from both these effects namely, random motions 
and coherent flows. In this paper we consider models for the redshift space distortions which combine both these effects. 
We compare the predictions of these models with the commonly used phenomenological model which has only the non-linear 
effects from random motions. We also compare all these models with N-body simulations and investigate which model best flts 
the N-body results. The different models are presented in Section 2 and the results of the comparison with N-body simulations 
are presented in Section 3. The galaxy pairwise velocity dispersion is a quantity which crops up in any discussion of the effects 
of redshift space distortions on the two-point correlation function. This quantity is very interesting in its own right and it 
has received a considerable amount of attention (Davis & Peebles 1983, Bean et al. 1983,Mo,Jing & Bornor 1993, Brainerd et 
al. 1994, SomerviUe,Primack & Nolthenius 1997,Bharadwaj 1997,Mo,Jing & Borner 1997, Landy, Szalay & Broadhurst 1998, 
Ratcliffe et al. 1998, Strauss,Ostriker & Cen 1998, Jing & Borner 1998,Jing & Borner 2001, Bharadwaj 2001, Sheth et al. 
2001, Del Popolo 2001 ). This quantity has been observationally determined for different redshift surveys(eg. Jing, Mo & 
Borner 1998, LCRS ; Zehavi et al. 2002, SDSS ; Hawkins et al. 2002 ,2dFGRS ). The models we use for the redshift space 
distortion also make definite predictions for the pair-wise velocity dispersion. In Section 4 we calculate the pair-wise velocity 
dispersion predicted by the different models and compare these with the pair-wise velocity dispersion determined directly 
from the N-body simulations. 

In Section 5 we discuss our results and present conclusions. 

We would also like to point out that the models which we have considered for ^' are very similar in spirit to those 
considered by Matsubara (1994), Regos & Szalay (1995) and Fisher (1995). 



2 MODELING 

The two-point statistics of the galaxy distribution in real space is completely quantified by the phase space distribution 
function P2(r, Vi,V2) which gives the probability density of finding a galaxy pair at a separation r, one member of the pair 
having peculiar velocity vi and the other V2. The redshift space two-point phase space distribution function p2(s, vi,V2) is 
related to its real space counterpart through 

/92(s, Vl,V2) = P2(s - nC/,Vi,V2) (9) 

where we have assumed the plane parallel approximation, and the units are chosen such that Hq = 1. Here U = n - (v2 — vi) 
is the line of sight component of the relative peculiar velocity of the galaxy pair. Integrating out the peculiar velocities gives 
us the redshift space two-point correlation function 
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l + r(s) = / p5(s,Vl,V2)d='vi d\2 (10) 



We next introduce the key assumption in the model, the assumption being that the peculiar velocity v of any galaxy can 
be written cis a sum of two parts v = vc + vh, where vc arises from large-scale coherent flows into overdense regions and out 
of underdense regions, and is a random part arising from galaxy motions in virialized clusters and other non-linear regions. 
The large-scale coherent flows are correlated with the density fluctuations which produce the flows, and the two are assumed 
to be related through linear theory. The two-point statistics of the coherent flow is quantified through the distribution function 
P2c(r, Vie, V2c) which is defined in exactly the same way as p2 the only difference being that p2c refers to only the part 
of the peculiar velocities which arises from the coherent flows. The statistical properties of the random part of the peculiar 
velocity are assumed to be isotropic and independent of the galaxy's location. Its' joint probability density can be written as 
P2ij(vii{, V2_r) = fl([vifl]a;)sr([vifl]y) . . . £r( [v2i?] z ) where [viij]^:, [vifl]y etc. refer to the different Cartesian components of vijj 
and V2i?, and g{vR) is the distribution function for a single component of the random part of a galaxy's peculiar velocity. The 
joint distribution of vi = vic +'Vir and V2 = V2C -|- V2h can be expressed in terms of the distribution functions for vic,V2C 
and vi_R, V2H as 

P2(r, Vl,V2)= J d^Vli{d^V2i{P2c(r, Vl - VlJ{,V2 - V2J{)p2H(viH,V2i{) (11) 

Using this in equations (9) and (10) to calculate the redshift space two-point correlation function we have 



1 -|- ^'^ (s) = / dui du2 



J dui du2 J 



d vid V2 P2c(s — n(?7 -|- W2 — Ml), vi, V2) 



g{ui)g(u2) (12) 



where Uj and U2 arc the lino of sight components of -vm and ■V2r respectively. The term in the square brackets [..] in equation 
(12) can, on comparison with equations (9) and (10), be identifled as the redshift space two-point correlation function if only 
the effects of the coherent flows are taken into account 

1 + Cc(s) = j d\ld\2P2c{s-hU,Vl,V2) (13) 

and can be expressed as 

^"(5)= J duidu2^h{s- 'n{u2 - ui)) g{ui) g{u2) . (14) 

To summarize, we start from the assumption that the galaxy peculiar velocities have two parts, one from the coherent flows 

and the other from random motions. We show that the redshift space correlation function is ^J, which has only the effect 
of the coherent flows, convolved along the line of sight with the one- dimensional distribution function of the random part of 
the galaxy's peculiar velocity, there being two convolutions, one for each galaxy in the pair. 

The fact that only the relative peculiar velocity v = U2 — uj between the two galaxies appears in equation (14) allows us 
to simplify it a little further. Equation (14) can be expressed it in terms of the self-convolution of g{vii) 

f{v)= J g{v - u) g{u) du . (15) 

The function f{v) may be interpreted as the distribution function for the line of sight component of the random part of the 
relative peculiar velocity v = U2 —u\ which is also called the pairwise velocity. Using this, we finally obtain (J' in terms of 
as 



r(s||,sx) = j dvec{sii+v,s±)f{v) (16) 

We now shift our attention to ^J, the redshift space two-point correlation function if only the coherent flows are taken 

into account. As mentioned earlier, we assume that we are working at large scales where linear theory holds in real space and 
the density fluctuations are a Gaussian random field. Expanding p2c(s — nil, vi, V2) in equation (13) in a Taylor series in the 
relative peculiar velocity U of the coherent flow we have 

l + ^c(s) = ^^i^ (^^) y"dVd't^2«7"p2c(s,vi,v2) (17) 



Retaining only the terms to order n = 2 we have 



Cl,(sx,sy) =C(s)- —Vp{s±,s^^) + --^ap{s±,sii) (18) 
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Table I. 

Model fi'") 

C ft. ^exp(^)(£ii+|.|) 

D ft. V3=^-P(i^) 



where ^ is the galaxy two-point correlation function in real space, Vp is the line of sight component of the mean relative 
velocity between the galaxy pair (also called mean pairwise velocity) 



r 



= --S||/?6(s) 



(19) 

(20) 



and (jp is the mean square of the line of sight component of the relative peculiar velocity (also called the pairwise velocity 

dispersion) 



o-p(sx,S||) = 



= 



d'^V\(fv2 U'^ P2C(S, Vl, V2) 



Us) 



(21) 



(22) 



Equation (18), combined with equations (20) and (22), is exactly the same as the linear redshift space two-point correlation 
function calculated by Hamilton (1992). Decomposing the angular dependence of equation (18) into Legendre polynomials one 
recovers exactly the same angular moments as equation (3), (4) and (5), and the odd moments and all even moments beyond 
I = 4 are zero. Using ^£ as given by equation (18) in equation (16) corresponds to the phenomenological model discussed 
earlier for the non- linear effects in and this is one of the models which we shall be considering in the paper. 

Going back to equation (17) for ^J, it is possible to exactly sum up the whole series keeping all powers of U (Bharadwaj 
2001). All the non- linear effects which arise due to the mapping from real space to redshift space are taken into account in 
this calculation, and the resulting redshift space two-point correlation function is given by 



1 + ^z,i(sx,S||) = y ds|| G(s||,o-p(sx,S|| -|-S||)) X 

S|lVp(sx,S|| -|-S|| 



(23) 



^r(sx,S|| +S||) + 1 ■ 



2a%{s±,s\^+s'u) I 4o-|,(sx,S|| -I-Sn) 



where we use 



G{x,a) 



1 



2'Ka 



■ cxp[- 



- — 1 



(24) 



to represent a normalized Gaussian distribution. 

We now have two different possibilities, ^£ or ^f,^, which we cam use for in equation (16) to calculate the full redshift 
space two-point correlation function . The function ^J. has only the effect of the coherent flows and it has to be convolved 
with f{v), the one dimensional distribution function for the random part of the pairwise velocity, to calculate In this paper 
we have tried out four different models which correspond to for different choices for and f{v). These are listed in Table I. 



To highlight the salient features of the four models. Model A uses ^f, for and an exponential for /(«). This is the 

phenomenological model discussed earlier. This model has been used extensively by different people when analyzing both 
N-body simulations and actual redshift surveys. Models B, C and D all use The difference between these models is in 
the choice of f{v). Model B uses an exponential form for f{v) and model D a Gaussian. Model C corresponds to a situation 
where the one dimensional distribution function for the random part of the galaxy peculiar velocity g{u) is assumed to be 
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an exponential. The function f[v) is now the convolution of two exponentials. All the models for f{v) have only one free 
parameter, cr^ which may be interpreted as the pairwise velocity dispersion of the random part of the peculiar velocity. 
In the next section we test the predictions of these models against the results of N-body simulations. 

3 RESULTS FOR f 

In this section we calculate ^" for the four models discussed earlier and compare the results against the predictions of N-body 
simulations. 

3.1 The N-body Simulations. 

We have used a Particle-Mesh (PM) N-body code to simulate the present distribution of dark matter in a comoving region 
[179.2 Mpc]^. The simulations were run using 256^ grid points at 0.7/i~^Mpc spacing with 128^ particles for a ACDM 
cosmological model with ftrno = 0.3, ^ao = 0.7 and h = 0.7. We have used a COBE normalized power spectrum with the 
shape parameter T = 0.2 for which as = 1.03. 

The low resolution N-body simulation used here is adequate for studying the deviations from the predictions of linear 
theory in redshift space on scales where the real space density fluctuations are well described by linear theory. We have 
restricted our analysis to scales larger than 5h~^ Mpc, though strictly speaking we would expect linear theory to be valid 
at scales larger than something like 8h~^ Mpc. To test that our low resolution simulations are not missing out any crucial 
feature either in real space or in redshift space, we have compared the results of our N-body simulations with the Virgo 
simulations (Jenkins et al. 1998) which have a higher resolution and a slightly different normalization with as = 0.9. We find 
that on the length-scales studied here, the results of our simulation arc consistent with the Virgo simulation both in real and 
redshift space. We show the results from the Virgo simulation alongside with those from our N-body simulation. Our N-body 
simulation was run for five independent realizations of the initial conditions. 

Assuming that galaxies trace mass, 10^ dark matter particles were chosen at random from the simulation volume and 
the entire analysis was carried out using these. The particle distribution in real space was taken over to redshift space in the 
plane parallel approximation. We determined the two-point correlation function for the particle distribution both in real and 
in redshift space. The angular dependence of the redshift space two-point correlation function was decomposed into Legendre 
polynomials, and the anisotropy in quantified through the ratios £.2{s)/^o{s) and C4(s)/Co(s)- We also estimated the ratio 
Q(s) (eq. 7) which is somewhat different from ^2(s)/^o(s) in the sense that it uses an integrated clustering measure instead 
of Co(s). This has the advantage that in the linear theory of redshift distortion the value of Q{s) is expected to be a constant 
irrespective of the shape of the real space correlation ^(s). Our simulations have 0^0 = 0.3 and 6=1 which corresponds to 
= 0.49, and we expect Q{s) = 0.57. 

The average and the la errorbars for ^o, Q ^4/^0 were calculated using the five realizations of our N-body 

simulations and the results arc shown in Figures 1 to 4 respectively. The points to note are 

(a.) The results of our simulation are consistent with those of the Virgo simulation which are also shown in the figures 
(b.) We see substantial deviations from the predictions of linear theory in redshift space on scales where it is known to hold 
in real space. This is best seen in the behaviour of Q{s) which is supposed to be a constant with value 0.57. We find that the 
value of Q is much below this even at scales as large as 20/i~^Mpc. The values of Q increases gradually toward the linear 
prediction all the way to length-scales as large as 30 — 40/i~^ Mpc where it finally appears to saturate at the linear prediction, 
(c.) The errorbars incrcas with increasing pair separation and they are quite large beyond 25 Mpc. We have tried using 
a larger number of particles to estimate but this does not reduce the errorbars leading to the conclusion that the we are 
limited by the cosmic-variance arising from the finite size of our simulation and not by Poisson noise. Larger simulations will 
be required to make more accurate predictions for the nature of the redshift space anisotropies. 

3.2 Fitting the models to N-body simulations 

All the models require the real space quantities ^{s), Vp{s±,s\\) and aj>{s±,s\\) as inputs to calculate S," in redshift space. 
We use ^, the real space correlation function averaged over five realizations of the N-body simulation, and its moments to 
calculate Vp{s±,s\\) and ap{sj_,s^^) using equations (20) and (22) respectively. Again, calculating using any of our models 
requires us to specify (3 and an. We have used (3 = 0.49 which is the value corresponding to the simulation parameters, and 
we treat an as a free parameter which we vary to obtain the best fit to the N-body results. For each model we fitted the 
model predictions for ^2(s)/^o(s) and Q{s) to the N-body results using a x'^ minimization with cr_R as the fitting parameter. 
There are good reason to believe that linear theory will not hold for s < 8h^^ Mpc and the fit was restricted to the region 
8 < s < 40 Mpc. To check if the models also work on length-scales which are mildly non-linear in real space, we have also 
carried out the fitting over the range 5 < s < 40 Mpc. 
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Table II. 



Model 


cr/{ (km/sec) 




<Ti{ (km/sec) 






S < ,s 10 h '-!\l]H- 




5 « < 101, -L\ii)r 




A 


684 


0.055 


760 


0.51 


B 


564 


0.054 


540 


0.19 


C 


520 


0.050 


452 


0.86 


D 


489 


0.114 


367 


2.20 



We find that for both ^2(s)/Co(s) and Q{s), the value of is minimized at nearly the same value of ctr, and so we quote 
the values only for Q{s). The best fit values of an and the corresponding values of Xmin per degree of freedom v are listed 
in Table II. The model predictions at the value of an which gives the best fit in the interval 8 < s < 40 Mpc are shown 
along with the results of our N-body simulations in figures 1 to 4. 

We find that all the models give a very good fit to the monopole (Figure 1), and the best fit predictions of the difi^crcut 
models are indistinguishable from one another. Considering next the anisotropies in (Figures 2 and 3) over the length-scales 
8<s< 40/1"^ Mpc, we find that all the models give a reasonably good fit. Model C has the smallest best fit x^/^) Models 
B, A and D follow in order of increasing 'x^ jv. The values are < 1 for all the models, indicating that all of them give acceptable 
fits. It should be noted that the best fit values of cr vary considerably across the different models, and Model A predicts a 
value considerably larger than the other models. Shifting our attention to the fits over the length-scales 5 < s < 40/i~^ Mpc 
we find that model B gives the lowest value of the best fit x^ /v, followed by Models A, C and D. All the models, except model 
D, have best fit /v below unity and hence give acceptable fits. Interestingly, the acceptable models A, B and C seem to 
work better than one would expect given the fact that the length-scales ~ 5 Mpc would be mildly non-linear in real space. 
Model D shows considerable deviations from the N-body results at length-scales 5 < s < Mpc. Here again, the best fit 

values of aR show considerable variations across the models. Also, for the same model, the best fit aR changes considerably 
when the fitting is done over length-scales 5 < s < 40ft~^ Mpc instead of 8 < s < 40/i~^ Mpc. This is particularly noticeable 
for Model D where best fit cr_R decreases by ^ 25% when the filing is extended to smaller length-scales. This change is ~ 10% 
for Models A and C, and ~ 5% for Model B. It should also be noted that for Model A, the best fit (jr increases when the 
fitting is extended to smaller length-scales, whereas the effect is opposite in all the other models. 

We now turn our attention to the hexadecapolc ratio ^4/^0 (Figure 4). Here again, for all the models we use the values 
of OR for which the model predictions for Q{s) give the best fit to the N-body results. The ratio £,aI£,o calculated with these 
values of ctr are shown in Figure 4. We find that in the range 10 — 22 /i~^Mpc the predictions of all the models fall below the 
N-body results, These deviations are within the Icr crrorbars and larger simulations are required before we can be really sure 
of the statistical significance of this effect. We have also tried fitting our models to the N-body results using a x^ minimization 
for ^4/^0 with aR as the free parameter. The best fit aR obtained this way are quite different from those obtained by fitting 
and Q(s) and we do not report these values here. This discrepancy may be indicating the inability of these models 
to adequately describe the hexadecapole {4, but further studies using larger N-body simulations with smaller errorbars are 
required to reach a definite conclusion. 



4 THE PAIRWISE VELOCITY DISPERSIONS 

The pairwise velocity dispersion is an important statistical quantity which sheds light on the clustering of matter in the 
universe. There are several approaches to determine the pairwise velocity dispersion on small scales from observations, for 
example, using the cosmic virial theorem (Peebles 1980, Suto & Suginohara 1991, Del Popolo 2001) or by modeling the 
distortions in the redshift-space correlation function (eg. Davis, Geller & Huchra 1983, Davis & Peebles 1983, Bean et al. 
1983, Mo,Jing & Borner 1993, Jing, Mo & Borner 1998, Jing & Borner 1998, Landy, Szalay & Broadhurst 1998, Ratcliffe et 
al. 1998, Zehavi et al. 2002, Hawkins et al. 2002). 

Our interest lies in the fact that the models which we have used to fit also make definite predictions for the pairwise 
velocity dispersion at large scales where we expect linear theory to hold. The pairwise velocity dispersion crfj, a symmetric 
rank two tensor, is defined as the second moment of the relative velocities of galaxy pairs and its value can be calculated from 
the distribution function p2(r, vi,V2) as 

'^%{r)= y(v2-vi)i(v2-vi)jP2(r,vi,V2)dVd^t;2/[l-h^(r)] (25) 

where i,j refer to different Cartesian components. Our work is restricted to large scales where linear theory holds in real space 
and we use 1 -|- ^(r) w 1. The most general form for a'^jijc) which is consistent with statistical homogeneity and isotropy is 
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Figure 1. This shows the monopolc as determined from our N-body simulations and the Virgo sinmlation. The normalization of the 
power spectrum used in the Virgo simulation is slightly different from the one used by us (Section 3.1), and the results from the Virgo 
simulation have been appropriately scaled to compensate for this. The figure also shows the predictions of the four models considered 
here for the value of ctjj (Table II) which gives the best fit to Q{s) in the interval 8 < s < 40 h~^Ts/lpc. The outcome of our simulations, 
the Virgo simulation and the best fit predictions of all four models are indistinguishable from one another. 
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Figure 2. This shows the ratio ^|/?o ^® determined from our N-body simulations and the Virgo simulation. It also shows the predictions 
of the four models considered here for the value of (Tr which gives the best fit in the interval 8 < s < 40 ?i~^Mpc (Table II). 



afjir) = (T'i{r)6ij + [<Tj|(r) - (T'i{r)]{rirj/r'^) . (26) 

Here cr'j^ir) is the pairwise velocity dispersion for the velocity component perpendicular to the pair separation r and o"y(r) 
is the dispersion for the velocity component parallel to r. The behaviour of crfj{r) is completely specified through these 
two components cr±{r) and o-jj(r). We next recollect the fundamental assumption underlying all the models which we have 
considered in the previous section ie. the peculiar velocity of any galaxy has two parts, one arising from coherent flows 
and another from random motions. Under this assumption the two-point distribution function p2 is the convolution of two 
distribution functions (eq. 11) one describing the two-point statistics of the coherent flow and another for the random motions. 
Using this in (equation 25) to calculate crfj{r) gives us 

<^l{r) = <Tlc{r) + 'r% (27) 

<rl{r) = (Tlc{r) + <T% (28) 

for all the models. Here cr^ is the isotropic contribution from random motions, and crfci''^) '''j.c(^) the contributions 
from coherent flows. 

Proceeding in exactly the same way as when using the models to fit S," , we assume that the coherent flows are related to 
the density fluctuations through linear theory ie. a]_c = a]_i^ and cr|p = cr|^. This allows us to express crj|^(r) and crxc(^) in 
terms of the moments of the real space two-point correlation function (Bharadwaj 2001) as 
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Figure 3. This shows Q{s) as determined from our N-body simulations and the Virgo simulations. It also shows the predictions of the 
four models considered here for the value of ctji which gives the best fit in the interval 8 < s < 40/i~^Mpc (Table II). The horizontal 
line at Q{s) = 0.57 is the constant value predicted by the linear theory of redshift distortions. 
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Figure 4. This shows the ratio ^|/i^q as determined from our N-body simulations and the Virgo simulations. It also shows the predictions 
of the four models considered here. It should be noted that model predictions are for the value of an which gives the best fit to Q{r) 
and not q/^g in the interval 8 < s < 40?i-iMpc (Table II). 



<rU^)=f3'r'[^i^{r)-^U{r)] (29) 

al^r) = 0' [^U'r) - ^6(r) + ^U'r)] ■ (30) 

In calculating ajj^^ and aj_j^ we have used the average real space two-point correlation function and its moments determined 
from our N-body simulations. In addition to crjj^ and crx^, all the models considered in this paper also need the value of an as 
an input to calculate and a^. In Section 3, for each model wo have determined the best fit value of an (Tabic II) for which 
the model predictions for Q{s) best match the N-body results in the range 8 < s < 40 Mpc. We have used these values of 
cTjj to calculate the pairwise velocity dispersion predicted by each of these models. The two independent components of the 
pairwisc velocity dispersion (cTj| and a±) were also determined directly from N-body simulations and the results are shown in 
Figure 5 and 6. 

We find that cry and crx determined from our N-body simulations decreases with increasing r at length-scales r < 
15 h~^Mpc, after which it is more or loss constant with possibly a very slow variation with r. Our N-body results are consistent 
with the high resolution simulations of Jenkins ot al. (1998). It is important to note that the variation of cr|j(r) with r plays 
an important role in redshift space distortions. For example, at linear order (eq. 18) the redshift space two-point correlation 
function $*(sx,S||) depends explicitly on ■^ap{s±,s\\) which is the second derivative of the line of sight component of the 

pairwise velocity dispersion. All the terms involving in the expressions for the different angular moments of (eqs. 4, 5 
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Figure 5. This shows (T|| as determined from our N-body simulation, along with the predictions of linear theory (eq. 29) and all the 
models considered in Sections 3. The models differ from the linear predictions in that they also have a contribution from random motions 
added in quadrature to the linear predictions (eq. 27). 
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Figure 6. This shows a± as determined from our N-body simulation, along with the predictions of linear theory (eq. 30) and all the 
models considered in Sections 3. The models differ from the linear predictions in that they also have an contribution from random 
motions added in quadrature to the linear predictions (eq. 28). 



and 6) arise from this. The spatial variation of afj{r) also plays an important role in determining in equation (23) where 
all the non- linear effects of the mapping from real to redshift space are taken into account. 

Turning our attention to the model predictions, wc first consider cr|| and a± calculated using only linear theory (eqs. 
29 and 30) with the real space correlations ^(r) and its moments determined from N-body simulations, as inputs. We find 
that these fall short of the values of a\\ and a± determined directly from N-body simulations. Also, the r dependence of a\\ 
and a± are quite different, with the N-body results decreasing and the linear predictions increasing with increasing r. At 
length-scales r > 25h^^ Mpc, the curves showing linear theory and the N-body results are approximately parallel, with the 
linear predictions being approximately 50km/s below the N-body results. 

The model predictions differ from the linear theory predictions in that they have a contribution from random motions ur 
added in quadrature to the linear predictions (eqs. 27 and 28). One might hope that the contribution from random motions 
will compensate for the shortfall in the linear predictions relative to the N-body simulations, and the predictions of the two 
will match at least at length-scales r > 25 Mpc whore the two curves are parallel. The problem is that all the models 
predict different values for an, and the predicted values are too large. Model A which has the highest value of an fares the 
worst with the predicted a\\ and a± being much larger than the N-body results. The predictions of Models B, C and D are 
slightly closer to the N-body results, but they arc all still very significantly higher than the N-body results. In summary (T|| 
and a± predicted by all the models are significantly in excess of the values determined directly from N-body siirmlatioris. This 
indicates that there is a gap in our understanding of what is really going on. 

The possibility of using the pairwise velocity dispersion as a tool for distinguishing between different cosmological models 
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has been controversial and this has been hotly debated in the literature (eg. Ostriker & Suto 1990, Cen & Ostriker 1992, 
Couchman & Carlberg 1992, Gelb & Bertschinger 1994, Zurek et al. 1994, Brainerd & Villumsen 1994, Brainerd et al. 1994, 
Somcrvillo,Primack & Nolthenius 1997). An important fact that wc should remember while measuring the pairwisc velocity 
dispersion from N-body simulations is that it is a pair weighted statistic and is heavily weighted by the densest regions 
present in the sample. These regions naturally have the highest velocity dispersion and this tends to push up the estimate. 
The statistic is strongly dependent on the presence or absence of rich clusters within the sample (eg. Mo,Jing & Borner 1993, 
Marzke et al. 1995, Mo,Jing & Borner 1997, Somerville,Primack & Nolthenius 1997, Guzzo et al. 1997, Hatton & Cole 1999). 
It has also been confirmed by several authors (eg. Sandage 1986, Brown & Peebles 1987, Willick et al. 1997, Strauss, Ostriker 
& Cen (1998)) that the velocity field is very cold outside the clusters. We note that these effects arc not very crucial in our 
work. This is because we have used exactly the same set of particles drawn from our N-body simulations to determine both ^" 
and afj , and we have been testing if the models which make reasonably good predictions for are also successful in correctly 
predicting tj^j . We would expect this to be true because the peculiar velocities which are quantified by the pairwisc velocity 
dispersion are also the cause of the redshift space distortions. Surprisingly, we find that the model predictions for a^j are 
significantly in excess of afj determined directly from the simulations. 



5 DISCUSSION AND CONCLUSIONS. 

The galax;y two-point correlation function determined from redshift surveys shows significant deviations from the predictions 
of the linear theory of redshift space distortion even on scales as large as 20 — 30 /i~^Mpc where linear theory is expected to be 
valid on real space. Any attempt to determine B from redshift surveys requires that these deviations be properly modeled and 
taken into account. Modeling redshift space distortions basically requires a joint model for galaxy peculiar velocities and their 
correlations with the galaxy clustering pattern. Such models test our understanding of the gravitational instability process by 
which the large scale structures are believed to have formed. 

We have considered four different models (details in Section 2) for the redshift space two-point correlation function ^" . All 
the models are based on the assumption that galaxy peculiar velocities may be decomposed into two parts, one arising from 
coherent flows and another from random motions. It is also assumed that in real space the coherent flows are well described by 
the linear theory of density perturbation. Deviations from the predictions of the linear theory of redshift space distortion arise 
from two distinct causes which aifect the mapping from real to redshift space (a.) non- linear effects due to the coherent flows 
(b.) the random motions. Among the four models, Model A does not incorporate the non-linear effects due to the coherent 
flows. It combines the predictions of the linear theory of redshift space distortion (Kaiser 1987, Hamilton 1992) with the effect 
of the random motions which is modeled through an exponential distribution function for the pairwise velocity. This is the 
popular phenomenological model which has been widely applied to the analysis of galaxy redshift surveys (eg. Hawkins et al. 
2002). Models B, C and D all take into account non-linear effects arising from the coherent flows (Bharadwaj 2001), and they 
differ from one another in the choice of the distribution function for the random part of the pairwise velocity. 

All the models have only one free parameter, an which is the one dimensional random pairwise velocity dispersion. For 
each model we have determined the value of an for which the model predictions best flt the quadrupole anisotropy of 
determined from N-body simulations. We find that Model C gives the lowest value of the best fit /'^ over the range of 
length-scales 8 < s < 40/i~^Mpc where we expect linear theory to be valid in real space. In this model the distribution 
function for the random part of a galaxy's peculiar velocity is modeled as an exponential function. It may be noted that the 
other three models also give acceptable fits to the N-body results. 

We find that three of the models (A, B and C) also give acceptable fits over length-scales 5 < s < 40fe~^Mpc which 
includes a small region where perturbations are expected to be mildly non-linear. Model D where the distribution function 
f{v) for the random part of pairwise velocity is a Gaussian fails to give an acceptable fit. Model B where f{v) is an exponential 
gives the lowest value of best fit 'x^ jv. It may be noted that though the best fit value of 'x^ jv for model A, the commonly used 
phenomenological model, is around three times larger than that for Model B, it is not possible to draw a strong statistical 
conclusion as to which model is superior. This is because x^/^" < 1 for Models A,B and C and they all provide acceptable 
fits. The present work is limited by the large statistical error-bars in the quantities determined at large scales from N-body 
simulations. These errors arise from the limited volume of the simulations (cosmic variance). It should be possible to achieve 
lower \—a error-bars using larger simulations whereby we could distinguish between these models at a higher level of statistical 
significance. Wc propose to carry this out in the future. 

Interestingly, the best fit value of or shows substantial variations across the models. The best fit value of or is substantially 
smaller in the models which incorporate the non-linear effects of the coherent flows (B, C and D) as compared to Model A 
which does not include these effects. This indicates that there are significant nonlinear effects in the mapping from real to 
redshift space arising from the coherent flows. The commonly used phenomenological model does not incorporate these effects 
and in this model all deviations from the linear predictions are attributed to random motions. This leads to the pairwise 
velocity dispersion of the random motion [or) being overestimated. For eaxample, Hawkins et al. (2002) have used Model A 
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to fit the redshift space two-point correlation function of tlie 2dFGRS to obtain tlie best fit value an = 506 ± 52km/s. The 
findings of this paper show that Models B and C would be equally successful in fitting the same observation, and the best fit 
value of aR would be different for each of these models. This raises questions as to the interpretation of ctr determined by 
this method. 

Although the models are all reasonably successful in fitting the quadrupole anisotropics of , the model predictions for 

the pairwisc velocity dispersion are much larger than the values determined directly from N-body simulations. Surprisingly, 
the predictions of linear theory which has a contribution from only the coherent flows and not the random motions are much 
closer to the N-body results as compared to the model predictions. At large scales the predictions of linear theory, all the 
models and the N-body results are all very similar. The linear theory predictions are slightly below the N-body results, and 
one would expect that it would be possible to recover the N-body results by combining the linear theory predictions with the 
contribution from random motions. Unfortunately, all the models appear to be overestimating the contribution from random 
motions and the model predictions are significantly in excess of the N-body results. Also, the predictions of Model A fare 
the worst in comparison to the other models. A possible explanation why equations (27) and (28) overpredicts the pairwise 
velocity dispersion is that the linear component of the peculiar velocity also makes a contribution to the random motions. It 
is possible that this is already present in aR, and it contributes more than its due share to the pairwise velocity dispersion. 

In the linear theory of redshift space distortions the hexadecapole anisotropy arises from the line of sight component of 
the pairwise velocity dispersion. The fact that none of the four models considered here give a very good fit to the hexadecapole 
is probably related to the fact that the models also do not predict the correct pairwise velocity dispersion. 

We note that the assumption that galaxy peculiar velocities can be decomposed into two parts, one coherent and another 
random is consistent with the halo model. The random part may be attributed to motions inside the halo and the coherent 
part to the overall motion of the halo. Seljak (2001), White (2001) and Kang et al. (2002) have calculated the galaxy power 
spectrum in redshift space using the halo model. Sheth & Diaferio (2001) have calculated the pairwise velocity dispersion 
using the halo model. 

In conclusion we note that the nonlinear effects in the mapping from real to redshift space introduced by the coherent 
flows are important. Models which incorporate these effects provide an equally good to the quadrupole anisotropies of as 
compared to models which are based on the linear theory of redshift distortion. Unfortunately, none of these models make 
correct predictions for the pairwise velocity dispersion. This indicates that there is a gap in our understanding of the statistical 
properties of the peculiar velocities and their effect on the redshift space two-point correlation function. 
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